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A search for — > v e oscillations has been conducted at the Los Alamos 
Meson Physics Facility using from ir + decay in flight. An excess in the 
number of beam-related events from the v e C — > e~ X inclusive reaction is 
observed. The excess is too large to be explained by normal v e contamination 
in the beam at a confidence level greater than 99%. If interpreted as an 
oscillation signal, the observed oscillation probability of (2.6±1.0±0.5) x 10 -3 
is consistent with the previously reported — ► v e oscillation evidence from 
LSND. 

14.60.Pq, 13.15.+g 
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I. INTRODUCTION 
A. Motivation 

In this paper we describe a search for neutrino oscillations from pion decay in flight 
(DIF). These data were obtained using the Liquid Scintillator Neutrino Detector (LSND) 
described in Ref. The result of a search for z/ M — > v e oscillations, using a flux from 
muon decay at rest (DAR), has already been reported in Ref. ||, where an excess of events 
was interpreted as evidence for neutrino oscillations. The present paper provides details of 
an analysis of the complementary process — > v e from neutrinos generated from tt + DIF. 

If indeed neutrino oscillations of the type z/ M — > v e do occur, then — > v e transitions 
must occur also. It is therefore important to search for the — > v e transition to demostrate 
that the DAR signal is due to oscillations, instead of being a property of the /i + decay. 
The 7r + DIF process provides a good setting for this search. It has completely different 
backgrounds and systematic errors from the DAR process, while providing an independent 
measurement of the same oscillation phenomena observed in the DAR measurement. Any 
excess of events in this analysis would support the neutrino oscillation hypothesis. 

The phenomenon of neutrino oscillations was first postulated by Pontecorvo || in 1957. 
The underlying theory has been described in detail in standard textbooks. A general formal- 
ism for neutrino oscillations would involve 6 parameters describing the mixing of all three 
generations and the possibility of CP violation. In general, a beam can oscillate into both 
v e and v T with different amplitudes and different distance scales, set by the three-generation 
mixing angles and the three mass-squared differences. In the present case a relatively pure 

beam is produced at the source. The LSND detector is sensitive to the v e state and 
thus, for simplicity, we approximate the process by a two-generation mixing model. The 
oscillation probability can then be written as 

P = sin 2 29 sin 2 ^1.27 Am 2 — ) , (1) 

where 9 is the mixing angle, Am 2 (eV 2 / c 4 ) is the difference of the squares of the masses of the 
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appropriate mass eigenstates, L (m) is the distance from neutrino production to detection, 
and E u (MeV) is the neutrino energy. The discussion is limited to this restricted formalism 
solely as a basis for experimental parameterization, and no judgement is made as to the 
simplicity of the actual situation. 

B. Comparison with other experiments 

In Ref. || the evidence restricting neutrino oscillation parameters is briefly reviewed. 
The salient features of that review are repeated here. There have been a series of experi- 
ments using beams derived from pion DIF which consist dominantly of with a small v e 
contamination. The most sensitive experiment was at Brookhaven in a specifically designed 
long baseline oscillation experiment, E776 [|]]. This limit is shown in Figure [I] along with the 
favored region obtained by the LSND experiment. The limiting systematic error in E776 
is a photon background from 7r° production, where one 7 is misidentified as an electron 
and the second 7 is not seen. The CCFR experiment || provides the most stringent limit 
on — > v e oscillations near Am 2 ~ 350 eV 2 /c 4 , but their limits are not as restrictive as 
E776 for values of Am 2 < 300 eV 2 /c 4 . The KARMEN experiment || has searched for 
v e oscillations using neutrinos from pion DAR. These neutrinos are monoenergetic, 
and the signature for oscillations is an electron energy peak at about 12 MeV. This method 
has very different backgrounds and systematics compared to the previous experiments but, 
unfortunately, does not yet have statistical precision sufficient to affect the exclusion region 
of Figure |I[ The KARMEN experiment also has searched for z/ M — > u e oscillations and has 
produced the exclusion plot shown in Figure |I[ This is currently the most sensitive limit 
experiment in this channel. KARMEN is located 18 m from the neutrino source, compared 
with 30 m for LSND. The experiments have sensitivities, therefore, that peak at different 
values of Am 2 . 

The most sensitive experiment searching for v e disappearance is Bugey J7| using a power 
reactor which is a prolific source of v e . The detectors at Bugey observe both the positron 
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from the primary neutrino interaction and the capture energy (4.8 MeV) from neutron 
absorption on 6 Li. The resulting limit is also shown in Figure |I[ 

The most sensitive searches for disappearance have been conducted by the CDHS 
|| and CCFR experiments. In each case two detectors are placed at different distances 
from the neutrino source, which is a DIF beam without focusing. The limits obtained by 
these experiments exclude the region with sin 2 29 > 0.08 for values of Am 2 typically above 1 
eV 2 /c 4 and are not as restrictive as the limits set by the appearance experiments described 
above. Finally, the E531 Fermilab experiment || searched for the appearance of tau decays 
from charged-current interactions in a high energy neutrino beam. This — > v T oscillation 
search excludes the region with sin 2 29 > 0.005 for values of Am 2 above approximately 10 
eV 2 /c 4 . Recently, the CHORUS and NOMAD experiments at CERN have reached limits 
close to that set by the E531 experiment with only a fraction of the data analyzed and 
should reach sensitivities of the order of 3 x 10 -4 in sin 2 29 in the near future. 

C. Experimental method 

LSND was designed to detect neutrinos originating in a proton target and beam stop 
at the Los Alamos Meson Physics Facility (LAMPF), and to search specifically for both 
Vfj, — > v e and — > u e transitions with high sensitivity. This paper focuses on the second of 
these two complementary searches. The neutrino source and detector are described in detail 
in Ref. ]lj, with a summary in Section II of this paper. For the DIF experimental strategy 
to be successful, the neutrino source must be dominated by u^, while producing relatively 
few u e by conventional means in the energy range of interest. The detector must be able to 
recognize v e interactions with precision and separate them from other backgrounds, many 
not related to the beam. The v e from conventional sources are small in number and are 
described in detail in Section VII. 

LSND detects u e via the inclusive charged-current reaction v e C — ► e~ X. The cross 
section for this process has been calculated in the continuum random phase approximation 
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(CRPA) fLOf . This calculation successfully predicts the v e C e~ X cross section from the 
fj+ DAR u e flux as measured by the LSND Q, KARMEN fl2j, and E225 |3| experiments. 
A similar calculation however predicts too large a cross section for the process v^C — > X 
at higher energies. A discussion of the cross section uncertainties and comparisons to the 
data is presented in Section VIII. The final state electron energy can range from zero to 
the incident neutrino energy minus 17.3 MeV, which corresponds to the binding energy 
difference between the initial nucleus and the final state nucleus in the ground state. 

The oscillation search analysis uses the following strategy. Beam-unrelated backgrounds 
induced by cosmic-ray interactions are removed as much as possible by requiring a positive 
identification of the electron from the v e C — > e~ X reaction in the tank. The remaining 
beam-unrelated background events in the sample are subtracted by using the data taken 
while the beam is off (beam-off sample) to determine the level of such background. Notice 
that the beam-off data is very well measured as LSND records approximately 13 times 
more data while the beam is off than while it is on. This procedure yields the number 
of excess events above cosmic background due to beam-induced neutrino processes. The 
remaining beam-related backgrounds are then subtracted to determine any excess above the 
expectation from conventional physics. The number and energy distribution of the excess 
events are used to determine a confidence region in the (sin 2 29, Am 2 ) parameter space. 

This paper describes two independent analyses that use independent reconstruction tech- 
niques and different event selections. This has allowed cross checks on the software and 
selection criteria and has resulted in a more efficient final event selection. They shall be 
referred to as "analysis A" and "analysis B" throughout this paper. 



D. Outline of the paper 

We present a brief description of the neutrino source and detector system in Section II. 
Section III describes the initial data selection for the DIF analysis, while the reconstruc- 
tion algorithm and particle identification parameters are discussed in Section IV. Section V 
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describes the event selection and efficiencies for two independent analyses. Distributions of 
the data are shown in Section VI. Section VII contains an assessment of the beam-induced 
neutrino backgrounds. Fits to the data and an interpretation of the data in terms of neutrino 
oscillations are presented in Section VIII. The conclusions are summarized in Section IX. 

II. NEUTRINO BEAM, DETECTOR AND DATA COLLECTION 

A. The neutrino source 

This experiment was carried out at LAMPF Q using 800 MeV protons from the linear 
accelerator. Pions were produced from 14772 Coulombs of proton beam at the primary 
beam stop over three years of operation between 1993 and 1995. There were 1787 Coulombs 
in 1993, 5904 Coulombs in 1994, and 7081 Coulombs in 1995. The fraction of the total DIF 
neutrino flux produced in each of the three years was 12% in 1993, 42% in 1994, and 46% 
in 1995. The flux in 1995 was slightly reduced with respect to the Coulomb fraction due to 
variations in the target conditions, which are described below. The duty ratio is defined to 
be the ratio of data collected with beam on to that with beam off. It averaged 0.070 for 
the entire data sample, and was 0.072, 0.078, and 0.060 for the years 1993, 1994, and 1995, 
respectively. 

A detailed description of the neutrino flux calculations from pion DIF in the LAMPF 
beam is given in Ref. fl4[]. A 1 mA beam of protons on the Al, A2 and A6 targets produces 
pions that are the source of the DIF neutrino beam |lj]. The primary source of neutrinos 
consists of a 30-cm long water target (A6) surrounded by steel shielding and followed by 
a copper beam dump. It is located approximately 30 m from the center of the detector. 
About 3.4% of the generated 7r + decay in flight due to the open space between the water 



The accelerator was operated under the name LAMPF until October 1995 when the name was 
changed to LANSCE (Los Alamos Neutron Scattering Center). 
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target and the beam stop, producing a flux with energies up to approximately 300 MeV. 
Most of the positive pions that decay come to rest prior to decaying. They then decay 
through the DAR sequence that produces the DAR neutrino fluxes via tt + — ► and 
fi + — > e + u e P l j / , where the v e and z/ M have a maximum energy of 52.8 MeV. Of the fi + that 
decay, approximately 0.001% decay in flight and produce a small v e contamination of the 

beam. Another small contamination comes from the decay mode ix + — > e + u e with a 
branching ratio of 1.24 x 10~ 4 . Together, these sources of v e constitute the major i/ e -induced 
background for the DIF analysis, as discussed in Section VII. The negative chain starting 
with tt~ leads to a smaller contamination of the beam with i? e , because the 7r~ production 
cross section is suppressed by a factor of about eight relative to the 7r + . 

The two upstream carbon targets Al and A2 were used to generate pion and muon 
beams for an experimental program in nuclear physics. They are located approximately 
135 m and 110 m, respectively, from the center of the detector. The flux from each target 
depended on the thickness as well as the proton energy in the primary beam reaching the 
target. They were originally 3 cm and 4 cm thick, respectively, and degraded slowly during 
the operation of the accelerator. Their thickness was monitored regularly during the runs 
and incorporated in the beam flux simulation. 

The DIF neutrino flux varies approximately as r~ 2 from the average neutrino production 
point, where r is the distance traveled by the neutrino. In addition, there is a significant 
angular dependence of the neutrino flux with respect to the direction of the incident proton 
beam. Thus, the DIF neutrino flux reaching the LSND apparatus has been calculated on 
a three-dimensional grid that covers uniformly the entire volume of the detector. The DIF 
fluxes at the detector center are illustrated in Figure |^ for the positive decay chains only. 
Figure 0(a) shows the flux from 7r + — > p^v^ while Figures 0(b) and 0(c) show the most 
significant v e background sources from 7r + — > e + u e and fi + — > e + u e u^, respectively. Notice 
that the contributions from the Al and A2 targets are generally small compared to that 
from A6. However, for z/ M — > v e oscillations with low Am 2 the flux from the two upstream 
targets can have a significant effect. 
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The systematic error on the DIF flux is estimated to be 15%. The calculated flux 
is confirmed within 15% statistical error by the LSND measurement of the exclusive 
Vpp — > pT 12 N gs reaction [TJ[] . This transition is very well understood theoretically, and 
the measurement is very clean due to the three-fold space-time correlations between the 
muon and the resulting decay electron and the positron emerging from the l2 N gs /5-decay. 
An independent beam flux simulation, based almost entirely on GEANT 3.21 [ 161 , has been 
developed in order to check the previous calculations and finds good agreement between 



calculated neutrino fluxes 17 



B. The detector and veto shield 

The detector consists of a steel tank filled with 167 metric tons of liquid scintillator 
and viewed by 1220 uniformly spaced 8" Hamamatsu photomultiplier tubes (PMT). The 
scintillator medium consists of mineral oil (CH 2 ) with a small admixture (0.031 g/1) of butyl- 
PBD. This mixture allows the detection of both Cerenkov and isotropic scintillation light, 
so that the on-line reconstruction software provides robust particle identification (PID) for 
e^, along with the event vertex and electron direction. The electronics and data acquisition 
(DAQ) systems were designed to detect related events separated in time. This is necessary 
both for neutrino-induced reactions and for cosmic-ray backgrounds. 

Despite 2.0 kg/cm 2 shielding above the detector tunnel, there remains a very large back- 
ground to the oscillation search due to cosmic rays, which is suppressed by about nine orders 
of magnitude to reach a sensitivity limited by the neutrino source itself. The 4 kHz cosmic- 
ray muon rate through the tank, of which about 10% stop and decay in the scintillator, is 
reduced by a veto shield to a 2 Hz rate. The veto shield encloses the detector on all sides 
except the bottom. Additional counters were placed below the veto shield after the 1993 run 
to reduce cosmic-ray background entering through the bottom support structure. The main 
veto shield [HJ consists of a 15-cm layer of liquid scintillator in an external tank, viewed 
by 292 uniformly spaced 5" EMI PMTs, and 15 cm of lead shot in an internal tank. This 
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combination of active and passive shielding tags cosmic-ray muons that stop in the lead shot. 
The veto shield threshold is set to 6 PMT hits. Above this value a veto signal holds off the 
trigger for 15.2 fis while inducing an 18% dead-time in the DAQ. A veto inefficiency < 10 -5 
is achieved off-line with this detector for incident charged particles. The veto inefficiency is 
larger for incident cosmic-ray neutrons. 

C. Detector simulation 

A GEANT 3.15-based Monte Carlo is employed to simulate interactions in the LSND 
tank and the response of the detector system [18| . It incorporates the important underlying 
physical processes such as energy loss by ionization, Bremsstrahlung, Compton scattering, 
pair production, and Cerenkov radiation. It also includes detector effects such as wavelength 
dependent light production, reflections, attenuation, pulse signal processing, and data ac- 
quisition. Much of the input to the detector response package was measured either in a 
test beam or in a controlled setting. Models for the transmission and absorption of light 
in the tank liquid were determined from measured data. The PMT characteristics, such as 
wavelength dependent quantum efficiencies, pulse shapes, and reflection characteristics, were 
measured and the results used in the simulation. The simulation is calibrated below 52.8 
MeV using Michel electrons from the decay of cosmic-ray muons that stop in the detector 
volume. The properties of the scintillator, including absorption length and detailed charac- 
teristics of Cerenkov radiation in this medium, are all checked in this way. The extrapolation 
to higher energies is then made using the MC simulation, which correctly incorporates the 
behavior of electrons in the detector medium. 

The primary Monte Carlo data set employed to calculate selection efficiencies is called 
the DIF-MC data set. The neutrino flux and energy spectrum were calculated at 25 points 
throughout the detector volume. The DIF-MC sample is created by folding the calculated 
Ufj, energy spectrum with the cross section predicted by the CRPA model to generate v e C —* 
e~ X interactions throughout the tank. This corresponds to 100% transmutation of the 
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beam to v e . The events were generated inside the surface formed by the faces of the PMTs. 

III. INITIAL DATA SELECTION 

The signature for the DIF oscillation search is the presence of an isolated, high-energy 
electron (60 < E e < 200 MeV) in the detector from the charged-current reaction u e C — > 
e~ X. The lower energy cut at 60 MeV is chosen to be above the Michel electron energy 
endpoint of 52.8 MeV, while the upper energy cut at 200 MeV is the point where the beam- 
off background starts to increase rapidly and the signal becomes negligible. The analysis 
relies solely on electron PID in an energy regime for which no control sample is available. 
Furthermore, with the exception of the v e C reaction leading to the 12 N ground state, there 
are no additional correlations that help improve the detection of the signal. 

The PID parameters used in the DAR analysis, Xa, Xr, Xt and Xtot - as defined in Ref. 
- have been used for the DIF analysis as an initial data selection. The disadvantages in this 
higher energy regime are that they do not discriminate adequately against a large beam- 
off background and are sensitive to energy extrapolation. Thus, the initial selection uses 
loose cuts based on the measured distributions of Xr, Xa and Xtot just below the Michel 
energy endpoint (50 < E e < 52 MeV), without any energy corrections. As demonstrated by 
MC simulations of the DIF data (DIF-MC), this selection is effective in identifying electron 
events from the v e C — > e~ X reaction. Over the energy interval of interest (60 < E e < 200 
MeV), the calculated efficiency is 98.1 ± 1.7%. 

In order to reduce the cosmic-ray muon induced background, we require the veto shield 
PMT hit multiplicity be < 4 for the data sample. In addition, the events were required to 
be reconstructed within a fiducial volume that extends up to 35 cm from the PMT faces. 
Space-time correlations have been used in the initial selection to reduce the background 
generated by the cosmic-ray muons, either directly or through the decay Michel electron. 
These correlations are described in the following subsections. 
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A. Future correlations 

Despite the veto shield hit multiplicity requirement, some cosmic-ray muons contaminate 
the sample. This is seen in Figure [3], which illustrates the distribution of the time difference 
between the current event and all the following ones, up to 51.2 fis. The fit to an exponential 
plus a constant reveals a time constant of 2.18 /is, identifying stopped cosmic-ray muons. 
Cosmic-ray muons that stop and decay in the detector are uniquely identified by the following 
Michel electron. As illustrated in Figure ^, there is a correlation between the muon tank hit 
multiplicity and the distance between the reconstructed vertices of the muon-electron pair. 
The difference in the samples shown in Figures |](a) and §(b) is briefly discussed below. 

Cosmic-ray muon events typically generate a high veto shield hit multiplicity. In order 
to suppress this high-rate background, as already mentioned in the previous section, the 
DAQ imposes a 15.2 /j,s dead-time after each event with a veto shield hit multiplicity > 6. 
Furthermore, all events with high veto shield hit multiplicities get a simpler event vertex 
reconstruction than the one described in Ref. [0, and no direction reconstruction. Correla- 
tions obtained from these data are shown in Figure |](a). The number of cosmic-ray muons 
with veto shield hit multiplicities < 6 is much lower than the ones with multiplicities > 6, 
which explains the smaller size of the sample shown in Figure §(b). Also, since these muons 
get both a full vertex and direction fit, the distance correlation between the muon-electron 
pair is tighter. In both cases, the distance correlation between the muon-electron pair de- 
grades with increasing muon tank hit multiplicity (or equivalently, energy) due to the on-line 
reconstruction algorithm which always assumes a point-like event. 

All events that are followed in the next 30 fis by an event with a tank hit multiplic- 
ity between 200 and 700 (typical for Michel electrons) are possible candidates for stopped 
cosmic-ray muons. If, in addition, the current event has a tank hit multiplicity above 600 or 
is reconstructed closer than 200 cm to the following Michel electron candidate, the current 
event is eliminated. The events that are eliminated by this selection are almost always fol- 
lowed by Michel electrons, as shown in Figure |^. This selection criterion is very powerful in 
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rejecting cosmic-ray muons and has a high efficiency for keeping candidate electron events 
from the v e C — > e~ X reaction. This efficiency is calculated to be 99.6 ± 0.4%. 

B. Past correlations 

Similarly, the time difference between the current event and all of the previous activities 
provides a distribution indicative of Michel electrons from stopped cosmic muons in the 
sample, as shown in Figure |^(a). Despite the energy requirement of at least 60 MeV there is 
still a small contamination from the tail of the Michel electron energy spectrum. Although 
this problem disappears at energies above 80 MeV, we choose to impose the following se- 
lection over the entire energy regime to maintain an energy independent selection efficiency. 
We require that the current event have no activities in the previous 30 /is with a tank hit 
multiplicity above 600 or closer than 200 cm. After imposing this cut, the time distribution 
with respect to previous events becomes flat, as shown in Figure ^|(b). Although this cut is 
powerful in rejecting the high-end tails of the Michel electron spectrum, this selection has 
an efficiency of only 85.5 ± 0.5%, due to the fact that it covers the entire energy interval. 

IV. EVENT RECONSTRUCTION AND ELECTRON IDENTIFICATION 

A. Introduction 

The event reconstruction and PID techniques that are used in the DIF analysis were 
developed to utilize fully the capabilities of the LSND apparatus. The basis for the recon- 
struction is a simple single track event model, parametrized by the track starting position 
and time (x,y,z,t), direction ((p,9), energy (E), and length (/). The coordinate system 
used throughout this analysis is located at the geometrical center of the detector, with the 
z-axis along the cylindrical axis of the tank (approximately parallel to and along the beam 
direction) and the y-axis vertical, pointing upwards. The expected PMT photon intensity 
and arrival time distributions for any given event are calculated from these parameters and 
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the result is compared with the measured values. A likelihood function that relates the 
measured PMT charge and time values to the calculated values is used to determine the 
best possible event parameters and at the same time provides PID. 

As mentioned in the introduction, two independent sets of reconstruction software were 
developed as a cross check of the analysis results. The two algorithms follow similar over- 
all strategies but differ in detail and implementation. The main differences lie in the 
parametrization of the various likelihoods and probability distributions that describe the 
detector response, and in the set of underlying event parameters used to describe the event. 

The electron identification is based on the relative likelihood of the measured PMT 
charges and times under the assumption that the source track is an electron. A detailed 
description of the physical processes in the tank can be found in Ref . [Q . Relativistic tracks 
in the detector generate light that falls into three categories: isotropic scintillation light that 
is directly proportional to the energy loss in the medium, direct Cerenkov light emitted in 
a 47° cone about the track direction, and isotropic scattered Cerenkov light. These three 
components occur in roughly equal proportions (0.35:0.32:0.33) for relativistic particles. 
Only the isotropic scintillation component occurs for non-relativistic charged particles. This 
difference forms the basis for distinguishing electrons from non-relativistic particles such as 
neutrons and protons. 

Each of the three light components has its own characteristic emission time distribution. 
The scintillation light has a small prompt peak plus a large tail which extends to hundreds 
of nanoseconds. The direct Cerenkov light is prompt and is measured with a resolution of 
approximately 1.5 ns. The scattered Cerenkov component has a time distribution between 
the direct Cerenkov light and the scintillation light, with a prompt peak and a tail that falls 
off more quickly than scintillation light. 

The two reconstruction algorithms used in the DIF analysis are based on maximizing 
the charge and time likelihood on an event-by-event basis. For any given event defined by 
the set of parameters a, 



14 



UCRHEP-E191/LA-UR-97-1998 

a = (x,y,z,t,(p,9,E,l), (2) 

the event likelihood for measuring the set of PMT charges (g,) and times (U) is written as 
a product over the 1220 individual tank PMTs as: 

1220 

Cevent = J[ Cq{qi\ «) Ct{U\ a) . (3) 
i=l 

Reversing the meaning of the likelihood function, C event is the likelihood that the event is 
characterized by the set a, given the set of measured charges (^) and times (£«). Maxi- 
mizing the event likelihood C event (or equivalently minimizing — \nC event ) with respect to a 
determines the optimal set of event parameters. 

The predicted likelihoods for PMT charges and photon arrival times are based on dis- 
tributions measured from a large sample of Michel electrons from stopped cosmic-ray muon 
decays, as described below. Analysis A uses the entire spectrum of Michel electrons, whereas 
analysis B uses only electrons with 38 < E e < 42 MeV, henceforth referred to as "monoen- 
ergetic". The upper edge of the Michel spectrum (52.8 MeV) is used to calibrate the energy 
scale of the system. The Michel electrons are well below the critical energy of 85 MeV and 
result in short track segments. The extension to longer, higher energy electron tracks is 
made by allowing for multiple discrete sources on the track. This is done either with two 
sources only along the track and fitting the distance between them (i.e., the track- length) in 
analysis A, or with a variable number of points, as determined from the energy of the event, 
distributed equidistantly along the track in analysis B. The energy dependence of the event 
likelihood is determined from the MC simulation. In the following subsections we briefly 
describe the charge and time likelihoods for a pointlike source of light. 

B. Charge likelihood 

Let us consider a pointlike light source located at (x, y, z) in the detector and direction 
determined by ((p, 9) in spherical coordinates. For low energy relativistic electrons the track 
length is comparable with the dimensions of the PMTs, and thus the pointlike approximation 
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provides a good model. The isotropic scintillation and scattered Cerenkov light have a 
combined strength <3> (photons per steradian), where $ is proportional to the energy E of 
the event. The strength of the anisotropic direct Cerenkov light is parametrized as p<3>, while 
the angular dependence is given by a nearly Gaussian function, /(cos0 e ). The angle 9 e is 
the angle with respect to the reconstructed event direction of the event and the function is 
normalized such that 

,,r /(cos0e)sin0 e d0 e = l. (4) 

The function /(cos0 e ), as determined from the data, is shown in Figure [F], with a vertical 
offset induced by the isotropic light. 

The average number of photoelectrons (PEs) /i expected at a phototube of quantum 
efficiency e, at a distance r from the source, and subtending a solid angle Q is given by 

fj, = e$F(cos0 e ,r) Q (5) 

in analysis A. The function F (cos 8 e ,r) is determined directly from the Michel data. In 
analysis B \x is parametrized as 

fi = e $ [e~ r/Xs n s + p /(cos 6 e ) e- r/Xc ft c ] . (6) 

The parameters A s and A c are the attenuation lengths for scintillation and direct Cerenkov 
light in the tank liquid, respectively. While these are indeed expected to be different for the 
different components of the light, the effective solid angles subtended by the PMT, Q s and fi c 
should, in principle, be identical. The difference is induced by the finite size of the PMTs and 



by the difference in the angular distributions of the two light sources |2(J • Both effective solid 
angles have been determined from the data. Notice that although the individual quantum 
efficiencies of the PMTs as well as the attenuation lengths are wave-length dependent, we 
use only global effective values as determined from the data. 

The probability of measuring n PEs in the presence of the light source is then given by 
a Poisson distribution of mean value /i, 
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P(n; fi) = i e-'V, (7) 

with /i given by Eq.([5]) or Eq.(|6]). However, since the LSND PMTs measure charge and not 
the number of PEs, the probability of measuring a charge q for a predicted value /i is given 
by 

oo 

V(q^) = Y, P ^n)P(n^), (8) 

n=0 

where the -P(g; n) functions are the charge response functions (CRFs) of the PMTs, i.e. the 
probability of measuring a charge q given a number of PEs n. Since 11 depends directly 
on the set of event parameters a, the probability V(q; fx) determines directly the charge 
likelihood £ q (q; a) for the PMT. 

In analysis A the V(q; n) functions are determined directly from the Michel data sample. 
In this sample, the predicted average number of PEs /i is calculated for every tube according 
to Eq.(]5p, for all of the events. For PMTs in a given predicted /x-bin, the distribution of the 
measured charge q, after proper normalization, yields directly the V(q; [/,) function required 
for the likelihood function above. The V(q; Li) functions obtained in this way contain all 
instrumental effects incurred in measuring the charge, such as saturation and threshold 
effects. Examples of these distributions are shown in Figure § for two predicted charges, \x 
= 0.0-0.5 PE and ll = 2.5-3.0 PE. 

Alternatively, analysis B obtains the two lowest CRFs, P(q; 0) and P(q; 1), and generates 
the higher P(q;n) distributions as follows. The lowest CRF, P(q;0), is just the Kronecker 
delta, 8 q> o, since the probability of measuring a charge q for no PEs vanishes identically for 
q > and is unity when q = 0. The second CRF, P(q; 1), is the single-PE response of the 
PMTs, as illustrated in Figure |9|. It is measured from low-intensity laser calibration data, 
taken during normal detector operation. The long tail of the single-PE charge distribution 
is probably due to collisions of electrons with material ahead of the first dynode. This effect 
is in good agreement with studies performed by the SNO experiment which uses the 



same type of PMTs. The higher CRFs (n > 2) are calculated by randomly sampling P(q; 1) 
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n times. With all CRFs normalized to unit area, the correct normalization of V(q; //) is 
automatically insured, 



Before going on to discussing the corrected time likelihood, we should point out that both 
reconstructions find a slightly longer attenuation length for the direct Cerenkov light in the 
MC Michel electrons than observed in the Michel electron data. This is needed to obtain 
agreement between the tank hit multiplicity and charge distributions obtained from the 
Michel electron data and the MC Michel electrons. This effect, which can be seen in Figure 
|7], has been shown not to affect significantly the charge and time likelihood distributions. 



The corrected time t c of a PMT is defined as the measured PMT time after corrections 
for the fitted event vertex time and light travel time from the event to the PMT surface. For 
prompt light this peaks at t c = with an RMS of approximately 1.5 ns. The time response 
functions for scintillation light and scattered Cerenkov light are more complicated and are 
determined from the Michel electron data. In addition, the time response functions depend 
upon the predicted charge. There is time slewing due to finite pulse rise time. There is also 
time jitter from the distribution of transit time of electrons in the PMT for signals with 
small numbers of PEs. Because the electronics responds to the first PE, above ten PEs the 
late tail in the distribution is negligible. Also, the amount of prompt Cerenkov light depends 
on whether or not the PMT is in the Cerenkov cone, as determined by f(cos8 e ). 

In analysis A the predicted average number of PEs \i is calculated for every PMT ac- 
cording to Eq.(|3p, for all events in the sample. For PMTs in a given predicted /i-bin, the 
distribution of the measured corrected time t c , after normalization to unit area, provides 
directly the probability P(t c ; /i) required for the time likelihood function, £ t (t;a). The 
P(t c ; n)s obtained in this way contain all instrumental effects incurred in measuring the 




(9) 



C. Corrected time likelihood 
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time, such as time slewing and PMT jitter. Figure [10] shows two examples of these distri- 
butions for Michel electron data and MC simulated data. Both distributions are obtained 
in the "isotropic" region cos# e < 0.3. Analysis A also measures the corrected time distribu- 
tions in the Cerenkov "peak" region (0.63 < cos# e < 0.73). An interpolation between these 
two distributions gives the time likelihood functions for the intermediate levels of direct 
Cerenkov light. 

The parametrization of the P(t c ; y) distributions of analysis B is described next. The 
timing distribution of the scintillation light for the LSND active medium has been measured 



to be of the form |22|] 



f(t)=A 1 e-*^ + 



A, 



for t > 0. 



(10) 



;i + t/r 2 ) 2 

The two terms above represent the fast and the slow components of the light, with time 
constants T\ = 1.65 ns and r 2 = 22.58 ns, respectively. The probability for observing a 
corrected time t r is 



P{t c 



N Jo 



f(t>) exp 



1 



A2 



2a 2 



dt', 



(11) 



which is the convolution of f(t) with a time smearing function, assumed to be a Gaussian of 
width o. The overall factor 1/N insures proper normalization to unity. While the integration 
of the fast component of the light can be analytically performed, this is not true for the slow 
component. Therefore we choose to parametrize the scintillation light as a superposition of 
three exponentially decaying functions, 



(12) 



Both the amplitudes and the time constants are also allowed to vary as a function of the 
predicted amount of scintillation light /i, as discussed above. The probability for recording 
a corrected time t r . is thus 

a 



N fr[ Jo [ 2cr 2 (n) 



1 3 

j7 22Mv) exp 

i=l 



2t 2 {h) nip) 



Erfc 



nil*) 
1 



dt' 

V^) 



(13) 
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with the normalization factor N given by 

iV = 2£AGu)T i (//). (14) 

i=l 

Replacing t c by t c — to(fi) in Eq.(|l^) above allows for additional time slewing corrections 
for the scintillation light. In LSND, the time slewing calibration is performed using laser 



calibration data p3fl , which is prompt. It is expected that the scintillation light should 
require additional time slewing corrections, which is confirmed by the data. A typical time 
probability distribution as a function of the corrected time is shown in Figure |ll|(a). The 
quality of the fit to the data is excellent and shows that the parametrization given by Eq.(|i~3"D 
is a very good approximation. 

The time distributions for the direct Cerenkov light are also measured from the "monoen- 
ergetic" Michel electrons for different values of the predicted charge fi. This is achieved by 
subtracting the appropriate underlying scintillation contributions from the corrected time 
distributions in the Cerenkov cone (0.53 < cos^ < 0.79). The resulting distribution is 
shown in Figure |TT|(b), which confirms the prompt character of the direct Cerenkov light. 

D. Electron identification 

The fitting procedures produce an accurate estimate of the amount of direct Cerenkov 
light in the event. In analysis A the level of Cerenkov light is determined after fitting the 
event to an electron model, i.e. a light source with an electron-like response for both charge 
and time likelihoods, that includes a full Cerenkov cone. With all other event parameters 
fixed, the amount of direct Cerenkov light is varied from none to the full amount for an 
electron event in order to maximize the event likelihood. This procedure determines a 
parameter Fc er which can thus vary between 0.0 (no direct Cerenkov light) and 1.0 (full 
amount of Cerenkov light). In analysis B the amount of Cerenkov light is determined by 
varying all event parameters including p, the Cerenkov-to-scintillation density ratio. Figure 
12] shows the distribution of the Cerenkov-to-scintillation density ratio for "monoenergetic" 
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Michel electrons, with a sharp peak at approximately p = 0.5. Moreover, for particles that 
are not expected to have any Cerenkov light (e.g. neutrons), the algorithm finds indeed a 
very low Cerenkov-to-scintillation fraction, which is still above zero because of fluctuations. 



This distribution is shown in Figure [13] for a sample of cosmic-ray neutrons in the same 
(electron-equivalent) energy range as the DIF sample. These have been tagged as neutrons 
by requiring the presence of a correlated 7 (from neutron capture on free protons, np — > dry) 
with a relatively high i? 7 parameter, as defined in Ref. f2j. Briefly, the _R 7 parameter is 
a quantity obtained from the 7 tank hit multiplicity, time and distance distributions with 
respect to the primary event. As shown in Ref. ||, it provides an excellent tool for identifying 
correlated photons and rejecting the accidental ones. 

The fitted optimum values of the charge and time negative log-likelihoods for the events 
are used as primary PID tools in the DIF analysis. In addition to the amount of Cerenkov 
light, they prove to be very different for non-electromagnetic events (e.g. neutrons) and 
provide very good discrimination against them, as will be shown in the next section. Analysis 
A makes use of the optimal values of the overall event charge and time likelihoods, C q and 
C t , and the likelihoods calculated in the Cerenkov cone only (0.53 < cos 6^ < 0.79), C qc and 
C tc . Analysis B uses instead the optimal values of the charge and time likelihoods obtained 
exclusively inside (LQ cer , Lxcer) and outside (LQ Sci , Lxsci) the cone. 

In both cases, the distribution of optimal likelihood values depend on the number of 
PMTs which were hit in a particular event. This is because the factor in the likelihood 
function for a PMT with a signal has a different functional form than for a tube without 
a signal. In order to remove this effect, the likelihoods are corrected as a function of the 
number of hit tubes. The mean value of the likelihood is then independent of the number of 
PMT hits. In addition, there is a dependence of the distribution of optimal likelihood values 



on the distance to the PMT wall, which is corrected for in an analogous way. Figure |H 
shows the corrected distributions for analysis A as obtained from the entire Michel electron 
spectrum. Distributions from analysis B, as obtained from "monoenergetic" Michel electron 
events, are illustrated in Figure [15|, before the hit multiplicity and distance corrections 
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described above. 

Both fitting algorithms significantly improve the position and direction accuracy over 



the angle resolution is approximately 6° for electron events over the energy interval of interest 
for this analysis. The energy resolution is limited to 6.6% at the Michel energy end-point, 
as stated in Ref . [U . This is due to the width of the single-PE response of the PMTs (Figure 
H) and also due to tube to tube variations in the response functions. 



The event selection presented in this section is designed to reduce cosmic-ray-induced 
background from the initial DIF sample. At the same time, the selection criteria must 
efficiently identify the electron of the final state in v e C — > e~X interactions. These events 
have the following characteristics: they have little or no activity in the veto shield system; 
they are nearly uniformly distributed inside the detector; they have no excess activity either 
before or after the event time; and they yield a track inside the tank which is consistent 
with an electron, as identified through the characteristic scintillation and Cerenkov light. 
These are the only features available for electron event selection except in the rare case of a 
transition to the 12 N ground state, which subsequently /3-decays with a 15.9 ms lifetime. 

The beam-off background data and simulated DIF-MC electron events are used in order 
to choose the optimal selection value for each quantity in a way unbiased by actual beam-on 
data. The sensitivity, or "merit" for the value of a selection parameter is defined as the 
efficiency e (determined from the DIF-MC) divided by the square root of the number of 
selected beam-off events, N ff, scaled by the duty ratio /: 



that used previously 0. The spatial position resolution is now approximately 11 cm and 



V. DATA SELECTION AND EFFICIENCIES 



A. Introduction 



M 



(15) 
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Each selection parameter value is varied and the point of maximum sensitivity, or "merit" , 
determines the optimal value of the selection. This method is independent of the beam-on 
data. Throughout this section the selection criteria for analyses A and B are discussed in 
parallel. 



B. Analyses A and B 

The cosmic-ray backgrounds are dominated by several types of processes. The level of 
all cosmic-ray-induced processes is measured in the beam-off data sample, and then the 
appropriate amount is subtracted from the final beam-on sample. 

Cosmic-ray neutrons generate a non-electromagnetic background. A fraction of these 
neutrons evade the veto shield and enter the detector volume to interact with carbon nuclei 
and protons in the liquid. The interaction length is approximately 75 cm. Their presence in 
the DIF sample is due to the very loose initial electron identification selection and can be 
consistently demonstrated in three different ways, as follows: 

The typical signature of neutron events is the 2.2 MeV correlated 7 that results from 
capture on free protons, np — > d'y. These 7 candidates are recorded in a 1000 /is interval 
after the primary trigger. The time difference between the primary events of the entire 



DIF sample (beam on + off) and all subsequent 7 events is shown in Figure |16|. The fitted 
lifetime of 186.2 ± 0.4 /is, is in very good agreement with the known neutron capture time 
of 186 fis. The constant part of the fit determines the total number of accidental photons. 
After subtracting this number from the total number of 7 candidates in the sample one 
obtains that on average every DIF event has one correlated photon. This is consistent with 
a considerable contamination of the DIF sample by cosmic-ray neutrons, since these are 
expected to have on average more than one associated 7. 



Secondly, as shown in Figure [H| neutron events do not produce significant Cerenkov light 
in the liquid. The distribution of the Cerenkov-to-scintillation density ratio, p, for the entire 
DIF sample (beam on + off) is shown in Figure IT7I for analysis B. The same distribution 
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for cosmic-ray neutron events with similar deposited energy and for the DIF-MC sample 
are superimposed. This illustrates that most of the (beam-off) background is dominated by 
non-electromagnetic events, consistent with neutrons. In order to select electron-like events, 
analysis A requires Fc er > 0.7 and analysis B requires p > 0.4, as dictated by the maximum 



merit algorithm. Figure [TJ] shows the Fcer and p distributions for the DIF sample and the 
DIF-MC electron events after all other selection criteria have been applied. 

Finally, the event charge and time likelihood parameters defined in Section IV are dif- 
ferent for neutron and electron events. These charge and time likelihood parameters are 
used in differentiating electromagnetic particles that produce Cerenkov light from non- 
electromagnetic backgrounds. Both analyses rely on this identification, using slightly differ- 
ent criteria. 

Analysis A uses a likelihood ratio, LR event , defined by forming the product of the charge 
and time likelihoods in each of the regions for the DIF beam-off sample and dividing it by 
the same product for the DIF-MC electrons: 

V(-\nC°/f) x V(-\nC°/J) x V(-\nC° t ff ) x V(-\nC°J f ) 
event ~ V(- In £™ c ) x V(- In £™ c ) x V{— In OT) x V(- In £? c c ) ' [ ' 

Figure ITU] shows the individual distributions of the four V(— lnC x ) functions for the beam- 



off data and for the DIF-MC electron data. The ratio LR event tends to be large for events 
that are like cosmic-ray background and small for electron-like events. Electron events are 
identified by requiring LR event < 0.5. The event time likelihood in the Cerenkov cone, L tc is 
also used in the selection, being sensitive to the presence of Cerenkov light. This parameter 
is required to have a value < 1.1. Analysis B uses only the individual time likelihoods for 
identifying electromagnetic particles, by requiring Lt sc % < 115 and Zr C er < 1-5. The charge 



and time negative log-likelihoods are illustrated in Figure ^0] for the entire DIF data (beam 
on + off), DIF-MC electron events and cosmic-ray neutron events. 

The second class of backgrounds is electromagnetic, and leads to events that are difficult 
to distinguish from pure electron events. Charged particles occasionally evade the veto shield 
and enter the liquid volume. These cannot travel into the liquid very far without depositing 
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large amounts of energy and are reconstructed with a position very close to the tank wall. 
Their reconstructed direction points predominantly into the detector and the track can be 
extrapolated back to the tank wall where the veto shield information can be used. The veto 
counter system that surrounds the detector provides PMT signals which are read out and 
recorded as are the tank PMTs. For events with a non-zero veto shield hit multiplicity, 
N ve to > 0, the reconstructed tracks in the detector are extrapolated backwards to intersect 
the veto shield. A corrected time difference, t veto , between the veto shield hit closest in time 
and space and the extrapolated event time is defined. Selecting events with \t veto \ > 50 ns 
(analysis A) or > 70 ns (analysis B) discriminates against any cosmic-ray-induced activity 
around the detector near the event in question. The t veto distribution for the entire DIF data 
sample (beam on+off) is shown in Figure In addition, a direct cut on the total veto hit 
multiplicity, N veto < 3, is required in analysis A. Analysis B does not impose this cut. Figure 
p2| shows the veto shield hit multiplicity distributions for the beam-on, beam-off and beam- 
excess events after all selection criteria of analysis B have been applied. The distribution 
for the beam-excess events is consistent with the veto shield accidental hit distribution. 

High energy 7 rays, from 7r° produced by neutron interactions in the lead shielding of 
the veto shield, enter the detector fiducial volume without leaving a veto signal. Energy is 
deposited through Compton scattering or by pair conversion. The latter process dominates 
above the 85 MeV critical energy of the liquid. The 7 attenuation length is roughly 50 cm, 
the radiation length in the liquid. The charged particles resulting from their interactions 
in the liquid point into the detector volume. These events are difficult to distinguish from 
electrons of the v e C — > e~ X reaction in this detector on the basis of electron identification 
alone. 

This class of backgrounds is characterized by its typical distribution of the length of 
the flight-path inside the detector. This quantity is defined as the distance between the 
reconstructed event vertex and the intersection of the backwards extrapolation along the 
reconstructed event direction with the PMT surface, A, (analysis A) or with the tank 
steel wall, S (analysis B). Although slightly different in their definitions, both quantities 
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correspond to the distance a neutral particle would have to travel in the liquid before it 
interacts. Events were required to satisfy A > 175 cm in analysis A and S > 225 cm in 



analysis B. The distributions of the A and S variables are shown in Figure |23] for the DIF 
sample and the DIF-MC electron events after all other selection criteria have been applied. 

Another possible source of cosmic-ray backgrounds that is reduced by the A and S 
selections arises from K° L decay inside the detector volume. The K° L can travel into the 
detector where the K e3 decay produces a positron (indistinguishable from an electron) and 
a 7r~ (K° L — > 7i~ e + v e ). The ir~ will stop and be absorbed, while the positron is in the energy 
range of interest. The other decay chain [K° L — > vr + e~ v e ) cannot contribute to the final DIF 
sample due to a degraded PID generated by the muon from the pion decay being virtually 
simultaneous with the electron, and due to space-time correlations with the positron from 
the subsequent muon decay. 

The electron from the v e C — > e~ X reaction is backward peaked, opposite the direction of 
the incident neutrino. Due to the geometry of the detector shielding, beam-off data favor the 
neutrino direction. Furthermore, the electron from one of the beam-induced backgrounds, 
the ve — > ve elastic scattering, is also strongly peaked in the forward direction. The cosine 
of the angle between the reconstructed event direction and the incident neutrino direction, 
cos 6^, is used to remove most of this background, by requiring cos^^ < 0.8 in both analyses. 

The veto system is very effective at rejecting cosmic-ray-induced backgrounds, but there 
were several penetrations in the system. A penetration at the lower upstream end of the 
veto system allows cables to enter the tank. For part of the data taking period there were 
several poorly performing PMTs at the top of the veto system. These regions were removed 
from the final data set of analysis A by requiring that the projected entry points of the 
events not lie in the regions (240° < 4> xz < 300°, -0.2 < cos9 y < -0.6) and (240° < <p xz < 
300°, 0.6 < cosOy < 1.0). The angles 6 y and (p xz are defined with respect to the coordinate 
system of the detector defined at the beginning of Section IV. 

The quantities A, F Cer , LR event , L tc , cos 6^, t veto , N veto and the "veto hot-spots" (analysis 
A) and S, p, Lxcer, Lxsci, cos 6 V and t veto (analysis B) are used to select a final sample of 
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events for the DIF analysis. The values of these selection criteria, efficiencies for z/ M — > 
v e events, and rejection power are shown in the Tables | and |I| for analyses A and B, 
respectively. Event selection efficiency is defined with respect to events generated inside the 
fiducial volume of the detector that extends all the way to the PMT surfaces id > cm). 

Table shows the number of beam-on, beam-off, and excess events that result from 
the event selections described above. There is a clear excess of events above beam-unrelated 
backgrounds that is consistent with v e C — > e~ X reactions. 

Average beam-unrelated backgrounds are determined by the total number of beam-off 
events times the beam duty ratio. A better estimate of the beam-unrelated background in 
the beam-on sample relies on further information contained in the beam-off event sample. 
The characteristic shape of the event direction distributions in the tank coordinate system 
for v e C — > e~ X signal events are different than those for the beam-off sample. Figure p4] 
shows the distributions in cos 9 y and (p xz of the event directions for the final beam-off DIF 
sample as compared to the same distributions as obtained for the DIF-MC electrons. 

It is possible to introduce a small bias into the event selection by using the maximum 
"merit" or sensitivity method to select values for the selection criteria. Because the algorithm 
uses the beam-off data, it can pick points where the beam-off data has fluctuated down. 
Even though this should be a negligible effect, the level of beam-unrelated backgrounds in 
the beam-on sample can also be determined nearly independently of the number in the beam- 
off sample. This is done by performing a maximum likelihood fit to obtain the number of 
beam-unrelated background events in the beam-on sample. The two-dimensional (cos 6 y , <fr xz ) 
distribution of the beam-on event directions is fitted to a sum of the shapes expected for 
signal events and beam-unrelated background events. The likelihood for the total number of 
beam-unrelated events is weighted by the Poisson probability expected from the predicted 
beam-unrelated average. The results of this procedure are shown in Table [TV] along with 
the results of using the product of the duty ratio and the number of beam-off events for 
each analysis. The probability that the number of observed beam-on events is a fluctuation 
is also shown in Table |TV|. A systematic uncertainty of 21% in the cross section, flux, and 
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efficiency is included in the calculation. 

The comparison of selections A and B along with the logical AND and logical OR of 
the two samples is shown in Table [V|. The number of beam-on events, background events, 
efficiencies, and resulting oscillation probabilities are all consistent within the statistical 
errors in the samples. Since the two analyses have low efficiencies, different reconstruction 
software, and different selection criteria, the overlap need not be large. The AND sample 
contains 8 beam-on events, which is consistent with the 11.6 events expected by comparing 
the overlap of DIF-MC data and beam-off data. The final event sample is obtained as 
the logical OR of the events from analysis A and analysis B. This procedure minimizes 
the sensitivity of the measurement to uncertainties in the efficiency calculations and also 
yields a larger efficiency than the individual analyses. The AND sample has the lowest 
efficiency for DIF electron events and therefore the least sensitivity. The OR sample has the 
largest efficiency and hence the largest sensitivity to oscillation signals. The probability that 
the backgrounds in the AND and OR samples fluctuate upward to the observed beam-on 
numbers are 0.18 and 1.1 x 10 -3 respectively. 



VI. DATA SIGNAL 



A. Distributions of data 



Extensive checks have been performed on the final DIF OR sample to study the con- 
sistency with electron events from the v e C — > e~ X reaction. The distribution of cos V , 
the cosine of the angle between the reconstructed electron direction and the incident neu- 



trino direction, is shown in Figure |25|. This distribution is slightly backwards peaked, as 
indeed expected from the u e C — > X reaction, and agrees well with that obtained from 
the DIF-MC. The distance to the PMT surfaces for the final beam-excess data set is shown 
in Figure |26|, which also agrees with the DIF-MC expected distribution. Notice that the 
apparent depletion of events in the outer region of the fiducial volume is caused primarily 
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by the A and S selections of the two analyses. Small deviations from the original (on-line) 
reconstructed event vertices, induced by the new reconstruction algorithms, contribute to 
a smaller extent to this effect. The x, y and z distributions for the final beam-excess DIF 
sample are shown in Figure [27J and are in very good agreement with those obtained from the 
DIF-MC simulations. The distributions of the 40 beam-on events and 175 beam-off events 



in the (x, y) and (y, z) planes are illustrated in Figure £8|. The energy distribution of the 
beam-excess events is illustrated in Figure ^9], together with the energy distribution of the 
beam-induced background and that expected from a positive — > v e oscillation signal for 
large values of Am 2 . 

B. Associated photons from neutron capture 

The v e C — > e~ X reaction is not normally expected to produce free neutrons at these 
energies. Only rarely (approximately 10%) is a neutron knocked out by the incident neutrino, 
and then it is identified by the presence of a correlated 2.2 MeV 7 from the capture on 
free protons. Also, the small u e contamination of the beam produces a small number of 
events with a correlated 7 via the inverse /3-decay reaction v e p — > e + n. The correlated 7 
identification relies on the _R 7 parameter mentioned earlier in the text, which in turn relies 
on the 7 tank hit multiplicity, time and distance distributions with respect to the primary 
event. The reconstruction algorithms used in the current analysis provide a better position 



resolution not only for the primary events, but also for the 7, as shown in Figure |30]. Using 
the sharper distance distribution between the 7 and the primary events in the calculation 
of Ry provides much better discrimination between correlated and accidental 7. 

The distributions for the number of photons with FLy > 1 are illustrated in Figure [31] 
for the final DIF beam-on, beam-off and beam-excess samples of analysis B. This particular 
value of the i? 7 cut accepts over 95% of the correlated photons, while at the same time 
rejecting approximately 95% of the accidental ones. The beam-induced excess yields a 
fraction of events with "correlated" photons (i? 7 > 1) that is consistent with that measured 
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in the v„ C 



H X channel, as reported in [15 



C. The transition to the 12 N gs 



The transition v e C — > e~ N gs , which is expected to occur roughly 5% of the time, 
is a useful signature in the search for — * v e oscillations. It is nearly free of cosmic-ray 
background due to the detection of the space-time correlated positron from the 12 N gs (3- 
decay yA \. It is noteworthy that the ground state of 12 N has been extensively studied in 



(p,n) reactions f25|, as well as its analog, 12 C (15.11 MeV), in (e,e') and (p,p') reactions. The 
transition is well known and can be characterized successfully. The positron has an end-point 
energy of 17.3 MeV and a decay time constant of 15.9 ms. The positron selection criteria 
are: (i) 0.052 ms < At < 45 ms; (ii) reconstructed distance to the primary electron < 100 
cm; (iii) tank hit multiplicity > 75 (in order to be above the accidental 7 ray background); 
(iv) positron energy < 18 MeV; (v) veto shield hit multiplicity < 4. Using the same selection 
criteria for the primary electron as described above, and in addition imposing the positron 
space-time correlations, 2 beam-on events and 1 beam-off event are observed. However, one 
of the two beam-on events has three correlated 7 and is thus not consistent with the ground 
state hypothesis. Eliminating this event from the sample, one obtains a beam-induced excess 
of 0.9 ± 1.0 events which is consistent with the expected 0.3-1.0 events, depending on the 
values of the oscillation parameters, as obtained from the inclusive analysis. 



VII. BEAM-RELATED BACKGROUNDS 

This section discusses the beam-related backgrounds induced by neutrino interactions, 
which are not accounted for by the beam-off subtraction. Each of the backgrounds depends 
upon neutrino fluxes and cross sections that are energy dependent. The proper variation of 
efficiency with energy in analyses A and B is used in Sections V and VIII. In this section 
a generic, energy independent efficiency for electron events of 10% is assumed for the sake 
of clarity, a number close to the actual average values for the two analyses. The effects of 
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systematic uncertainties on the beam-related backgrounds are crucial to the analyses. The 
major effects are discussed extensively in Section VIII. 

There are four significant neutrino backgrounds in the — > v e DIF oscillation search 
considered here. These backgrounds are: fi + — > e + v e v^ and 7r + — > e + v e DIF followed 
by v e C — > e~X scattering; ir + — > DIF followed by u^e —>■ v^e elastic scattering; 

and 7r + — > /i + z^ DIF followed by v^C — > v^Cn coherent scattering. Backgrounds from 
z/^C — > fi~X reactions are negligible. Muons that stop in the tank either decay or capture 
on carbon nuclei. The correlations in time and position between the muon and the secondary 
event removes them. In the rare case that the second event is missed, the lone muon fails the 
electron identification, and the event is almost always below the 60 MeV electron equivalent 
energy limit. In the case of fi~ decay in flight, the long lifetime of the muon and the electron 
identification requirements reduce this background to a negligible level. 



The significant backgrounds are summarized in Table [Vj for reconstructed event energies 
between 60 MeV and 200 MeV. The volume used for normalization throughout this section 
is the d > fiducial volume, as described earlier, which contains the equivalent of 5.4 x 10 30 
CH2 molecules (or 4.3 x 10 31 electrons). For the combined 1993, 1994, and 1995 running 
periods there were 9.2 x 10 22 protons on target (POT). 

The largest background is due to v e that come from /i + DIF in the beam stop, followed by 
v e C — > e~X scattering in the detector. This cross section is calculated in the CRPA model 
|T0| , as already mentioned. This results in a background of 3.8 events for the assumed 
efficiency The next largest is background is due to 7r + — > e + u e DIF in the beam stop, 
followed by u e C — > e~X scattering in the detector. The estimated background contribution 
is 1.9 events. The systematic error on this contribution is discussed in Section VIII. The 
previous two background catagories are produced by the same reaction as the DIF oscillation 
signal and are nearly impossible to distinguish from them on an event-by-event basis. 

Another background is coherent tt° production via the reaction v^C — > u^Cti . This 
cross section has been calculated in Ref. |2?J. Energetic electrons can be produced by the 
photons from the 7r° decay, which convert in the tank liquid and fake an electron. The 
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fraction of tt° that satisfy all selection criteria and are misidentified as electrons is 0.6. The 
estimated background contribution from coherent tt° production is 0.3 events. Note that 
the non-coherent ir° production is negligible at these energies. 

The last background considered is v^e —>■ v^e elastic scattering on the 4.3 x 10 31 electrons 
in the d > fiducial volume. This purely leptonic cross section is well known theoretically. 
The reaction is identified by the electron direction being nearly parallel to the incident 
neutrino direction. The fraction of electrons within the event selection region of cos6 u < 0.8 
is 0.05. The background contribution from v^e — > v^e elastic scattering is estimated to be 
0.1 events. 



Table [VT| shows the background estimates, where the total background is calculated to 
be 5.8 ± 1.0 events. The number of events expected for 100% — > u e transmutation is 4470 
events. 



VIII. INTERPRETATION OF THE DATA 

This section describes the interpretation of the observed event excess in terms of the 
theoretically expected processes and a neutrino oscillation model. The oscillation model 
employed here assumes two-generation mixing, as discussed in Section I. The confidence 
regions in the (sin 2 29, Am 2 ) parameter plane are calculated in this context. The effects of 
systematic errors are critical to the interpretation, and are described next. 



A. Systematic uncertainties 

The systematic uncertainties in this measurement arise from several sources. The dom- 
inant uncertainty comes from the knowledge of the underlying neutrino cross sections and 
the neutrino flux through the detector. The electron selection efficiency calculation also 
introduces some uncertainty. 

The oscillation search relies on the knowledge of the v e C — > e~ X cross section in the 
60-200 MeV electron energy range. The inclusive reaction has been calculated in the CRPA 
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model |1(J and has been measured [ 13(1 by using the v e flux from n + DAR. The DAR 



flux is in turn measured by the well understood ground state v e C — > e~ 12 N gs inverse f3- 
decay reaction. Both the ground state and the inclusive v e C reaction measurements agree 
well with the theoretical predictions, which indicates that both the flux and the cross section 
are predicted well. An extrapolation of the theoretical cross section must be made into the 
DIF energy region. We assign a systematic error of 10% to this extrapolation, based upon 
the agreement between the measured v e C data and the theoretical CRPA prediction. 

The DAR v e flux endpoint is at 52.8 MeV, below the region of interest for the DIF 
oscillation search. The DIF neutrino flux comes from pions that decay in flight rather than 
from stopped fi + . The u^C ground state and inclusive measurements of LSND |15] provide 
a check on the DIF flux. The v^C — > fi~ 12 N gs ground state cross section is also well 
understood and is nearly independent of energy above its 123 MeV threshold. Thus, LSND 
measures the integral of the flux above threshold. The agreement between the predicted 
flux and the measured flux gives a constraint on the flux above threshold with an error of 
15%. The C inclusive reaction cross section has a much stronger energy dependence than 
the ground state reaction. LSND measured this cross section with high statistics and 
obtained a value that is approximately 45% lower than the theoretically predicted value. 
The calculated flux x cross section does not agree with the measured data in this case. 
These neutrinos are in an energy range that overlaps with the DIF —>■ v e energy range 
and represent the same flux that the DIF oscillation search uses. It is possible that the 
flux x cross section for the u e C —>■ e~ X reaction also follows this trend and is lower than 
what we have assumed. The consequences of this are discussed below. 

The next important systematic error is the extrapolation of the electron identification 
efficiencies to energies above the Michel endpoint of 52.8 MeV. A GEANT 3.15-based Monte 
Carlo calculation is used for this purpose, as described in Section II. The MC generated 
events were checked against Michel data taken during the 1994-1995 run periods. The 
electron ID efficiencies are determined in the MC Michel electron sample and in the data 
Michel electron sample. The differences observed in these two samples result in a 15% 
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uncertainty in the selection efficiency. When the Michel MC sample is compared to the 
DIF-MC sample a lower difference of 12% is expected due to slightly narrower distributions 
in the DIF-MC sample. 

The effect of the systematic uncertainties on the oscillation search can be explained as 
follows. The DIF oscillation search looks for an excess signal in the v e C — > e~ X process 
above the background from the v e contamination in the beam. This background flux is pro- 
duced by the same DIF beam that produces the beam. The parent particle is dominantly 
either the /z + or the n + . The expected average number of events is given by 

Hiatal = e a Ue c {^^Pv^v e + $„J + Hbub (17) 

where e is the event selection efficiency, ^ Ue /u^ are the neutrino fluxes, a Ve c is the neutrino 
cross section, and [J, bub is the cosmic-ray background. The oscillation signal is proportional 
to the same product, £cr Ve c& Vli , as the neutrino background, since $ Ue is proportional to 
$ v . The effect of lowering the product e$a Ue c is to reduce the predicted beam-related 
background, i.e. the background from neutrino interactions from the v e contamination in 
the beam. This raises the observed oscillation signal. Only by raising the product e $ a Ve c 
is the oscillation signal decreased. 

In order to calculate conservative confidence regions, a value of the e$a Ve c 21% above 
the calculated value is assumed because of systematic uncertainties. This gives a larger 
confidence region because fewer of the excess events are attributed to oscillations. Then a 
value of e $ a UeC 45% lower than the calculated value is assumed, consistent with the LSND 
Vfj, C measurement. This has the effect of moving the lower contour to the right. The upper 
contour also moves somewhat to the right. The final confidence region is shown as the logical 
OR of the two extreme cases just described. 

B. Confidence regions 



In order to determine the significance of the observed signal in terms of potential neutrino 
oscillation effects, a confidence level calculation is made in the context of a two-generation 

34 



UCRHEP-E191/LA-UR-97-1998 

neutrino mixing model, as discussed in the Introduction. The oscillation probability is a 
function of the neutrino energy and the distance to the neutrino source. In the present case 
the distance to the source is ambiguous because of the presence of multiple beam targets, 
Al, A2, and A6. Therefore, the energy distribution alone is used to determine the confidence 
levels in the (sin 2 29, Am 2 ) parameter space. 

The data are binned into four equal energy bins between 60 MeV and 200 MeV. In each 
bin the DIF-MC data are used to calculate the expected number of oscillation events, fi osc , 
and beam-related background events, i^brb, at each (sin 2 29, Am 2 ) point. This number is 
added to the expected beam-unrelated background, Hbub, to determine the total expected 
number of events in each of the four energy bins: 

/7(sin 2 29, Am 2 ) = /7 osc (sin 2 29, Am 2 ) + fl BRB + jx B uB- (18) 

— * 

From the expected numbers, a four-dimensional Poisson probability density, p(N; fx) (one 
dimension for each bin) of all possible results for this experiment is determined. An inte- 
gration over all points in this space with a probability density greater than or equal to the 
measured data point value, 

V(N;fi) > V(N meas ;fI), (19) 
gives a probability for the (sin 2 29, Am 2 ) point: 

p(N;ft) where N= (i, j,k,l), i,j,k,l = 0,...,oo. (20) 

p{N \P)>p{N meas ]fl) 

This calculation determines confidence regions, or contours of equal probability, in the 
(sin 2 29, Am 2 ) space. As discussed above in Subsection A, the calculation is made for two 
extreme cases of the product e a Ve c The contours that result from the logical OR of these 
extremes are shown in Figure |32[ The calculation shows that the DIF result of this paper is 
consistent with the previous LSND DAR result |§ in terms of the two-generation oscillation 
parameters. 
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IX. CONCLUSIONS 

This paper reports a search for v e C — > e~ X interactions for electron energies 60 < E e < 
200 MeV. Table [VI] lists the expected contributions from conventional sources. This search 
is motivated by a high sensitivity to neutrino oscillations of the type — > v e , due to the 
small contribution from conventional processes to the u e flux in this energy regime. Two 
independent analyses observe a number of beam-on events significantly above the expected 
number from the sum of conventional beam- related processes and cosmic-ray (beam-off) 
events. The probability that the 12.5 (14.9) estimated background events fluctuate into 23 
(25) observed events is 7.0 x 10~ 3 (1.6 x 10~ 2 ). The excess events are consistent with — > u e 
oscillations with an oscillation probability of (2.6 ± 1.0 ± 0.5) x 10~ 3 . A fit to the event 
distributions, assuming neutrino oscillations as the source of v e , yields the allowed region 
in the (sin 2 29, Am 2 ) parameter space shown in Figure |32|. This allowed region is consistent 
with the allowed region from the DAR search reported earlier. This —>■ v e DIF oscillation 
search has completely different backgrounds and systematic errors from the —* v e DAR 
oscillation search and provides additional evidence that both effects are due to neutrino 
oscillations. 
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TABLE I. Selection criteria for analysis A. For each criterion is listed: the value of the criterion; 
the efficiency of the criterion after all other criteria have been applied; the number of events rejected 
by that criterion after all other criteria have been applied; the number of beam-on, beam-off and 
beam-excess events prior to applying the criterion. 



Criterion 


Cut Value 


Efficiency 


N ■ 


On 


Off 


Excess 


Fcer 


> 0.7 


0.68 


643 


60 


720 


9.6 ± 7.7 


A 


> 175 cm 


0.56 


516 


53 


600 


11.4 ± 7.3 


LRevent 


< 0.5 


0.89 


118 


32 


223 


16.7 ± 5.7 


Veto time 


> 50 ns 


0.98 


27 


26 


138 


16.4 ± 5.1 


cos6 u 


< 0.8 


0.95 


22 


24 


135 


14.7 ± 4.9 


Hot Spots 


(see text) 


0.94 


22 


25 


134 


15.7 ± 4.9 


Ltc 


< 1.1 


0.98 


11 


23 


125 


14.4 ± 4.8 


veto hits 


< 3 


0.98 


3 


23 


117 


14.9 ± 4.8 
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TABLE II. Selection criteria for analysis B. For each criterion is listed: the value of the 
criterion; the efficiency of the criterion after all other criteria have been applied; the number of 
events rejected by that criterion after all other criteria have been applied; the number of beam-on, 
beam-off and beam-excess events prior to applying the criterion. 



Criterion 


Cut Value 


Efficiency 


N re j 


On 


Off 


Excess 


S 


> 225 cm 


0.47 


1009 


89 


1037 


16.8 ± 9.7 


P 


> 0.4 


0.83 


689 


68 


738 


12.7 ± 8.5 


LTsci 


< 1.15 


0.90 


131 


34 


214 


19.3 ± 5.9 


Veto time 


> 70 ns 


0.96 


36 


27 


126 


18.0 ± 5.3 


cos6 u 


< 0.8 


0.95 


14 


27 


104 


19.7 ± 5.2 


LTcer 


< 1.5 


0.98 


10 


26 


101 


18.9 ± 4.1 
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TABLE III. Event count after all selection criteria have been applied (analyses A/B). Qag is 
the number of protons on target in Coulombs. 



Excess (A/B) 



Year 


QA6 


Beam On (A/B) 


Beam Off (A/B) 


Duty ratio 


1993 


1787 


1/2 


17 / 21 


0.072 


1994 


5904 


12 / 15 


42 / 41 


0.078 


1995 


7081 


10/8 


55 / 30 


0.060 


Total 


14772 


23 / 25 


114 / 92 


0.070 



-0.2 ± 1.0 / 0.5 ± 1.5 
8.7 ± 3.5 / 11.8 ± 3.9 
6.7 ± 3.2 / 6.2 ± 2.8 
15.1 ± 4.9 / 18.5 ± 5.0 



TABLE IV. Backgrounds and observed numbers of beam-on events. The beam-unrelated 
cosmic backgrounds are calculated in two ways. The product of the beam duty ratio and the number 
of beam-off events is shown first. The correponding result for the fitted number of beam- unrelated 
background events in the beam-on sample is shown in parenthesis. The distinction is described in 
the text. 



Analysis A Analysis B 



Beam-unrelated background 8.0 ± 0.7 (6.2 ± 2.0) 6.4 ± 0.7 (5.9 ± 1.9) 

Expected beam-related background 4.5 ± 0.9 8.5 ± 1.7 

Total expected background 12.5 ± 1.1 (10.7 ± 2.2) 14.9 ± 1.8 (14.4 ± 2.6) 
Observed beam-on events 23 25 

Fluctuation probability 7.0 x 10 -3 (1.2 x 10~ 3 ) 1.6 x 10~ 2 (1.2 x 10~ 2 ) 
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TABLE V. Comparison of results for the A, B, AND, and OR data sets. The BUB is the 
beam-unrelated background from cosmic rays. 



Sample 


Beam On/Off 


BUB 


^-Background 


Osc. Excess 


Efficiency 


Osc. Probability 


Analysis A 


23/114 


8.0 ±0.7 


4.5 ±0.9 


10.5 ±4.9 


0.084 


(2.9 ± 1.4) x 10~ 3 


Analysis B 


25/ 92 


6.4 ±0.7 


8.5 ± 1.7 


10.1 ±5.3 


0.138 


(1.7 ±0.9) x 10~ 3 


AND 


8/ 31 


2.2 ±0.3 


3.1 ±0.6 


2.7 ±2.9 


0.055 


(1.1 ± 1.2) x lO- 3 


OR 


40/175 


12.3 ±0.9 


9.6 ± 1.9 


18.1 ±6.6 


0.165 


(2.6 ± 1.0) x 10~ 3 



TABLE VI. The background estimates for the — > v e oscillation search are shown for a d > 
detector fiducial volume, 9.2 x 10 22 protons on target, and for reconstructed energies between 60 
MeV and 200 MeV. The number of events for 100% v„ — > v e transmutation is shown also. These 
numbers are illustrative for an electron selection efficiency of 0.10, independent of energy. The 
actual efficiencies in analyses A and B are slightly different and energy dependent. 



Process 


Flux (cm- 2 /POT) 


<a> u (10~ 40 cm 2 ) 


Eff. 


Number of Events 


v e C -► e~X {n DIF) 


3.8 x 10~ 14 


28.3 


0.10 


3.8 


v e C -► e~X (?r DIF) 


8.3 x 10~ 15 


79.2 


0.10 


1.6 




6.5 x lO" 11 


1.6 


0.06 


0.3 


v^e -> v^e 


6.5 x lO -11 


0.00136 


0.005 


0.1 


Total background 








5.8 


100% Ufj, — > v e transmutation 








4470 
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FIGURES 




FIG. 1. The LSND (sin 2 29, Am 2 ) favored regions obtained from the — > v e DAR oscilla- 
tions search. The darkly-shaded and lightly-shaded regions correspond to 90% and 99% likelihood 
regions. Also shown are the 90% confidence level limits from KARMEN (dashed), E776 (dotted) 
and the Bugey reactor experiment (dot-dashed). 
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FIG. 2. Calculated and v e DIF fluxes at the detector center from the A6 target (solid 
histograms) and from the A1+A2 targets (dashed histograms). 
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FIG. 3. Time difference distribution to events subsequent to all primary events for the initial 
DIF data after standard PID, veto shield hit multiplicity and fiducial volume cuts. 
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FIG. 4. Reconstructed distance distribution between Michel electron and parent stopped 
cosmic-ray muon, versus muon tank hit multiplicity. Figure (a) is for muons with a veto shield 
hit multiplicity > 6 and (b) is for muons with a veto shield hit multiplicity < 6. The regions in 
the upper left corners, delimited by the dashed lines, are the regions allowed by this selection (see 
text) . 
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Future event dt ((is) Future event E (MeV) 

FIG. 5. (a) Time difference to future events and (b) energy distribution of the future events 
for primary events of the initial DIF sample that fail the future space-time correlations. The fit 
in (a) is to an exponential plus a constant with a time constant of 2.2 fj,s. The fit in (b) is to the 
Michel electron spectrum shape. 
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FIG. 6. Time difference to all previous events (a) before and (b) after the past space-time 
correlation cuts. The sharp edge at 15.2 /xs in both (a) and (b) is a reflection of the DAQ operation, 
as described in the text. The sharp edge at 30 /is in (b) is induced by the selection algorithm. 
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FIG. 7. The ratio of the average measured charge to the predicted scintillation charge versus 
cos 9 e , the cosine of the angle between the event direction and the PMTs. The solid histogram 
shows the data, and the points with error bars show the MC simulation. 
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FIG. 8. Unnormalized charge response functions of the PMTs for two values of the predicted 
charge (a) /x = 0.0-0.5 PE and (b) \i = 2.5-3.0 PE. The solid histograms show the data, and the 
dashed histograms show the MC simulation. 
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FIG. 9. Single-PE charge response function of the PMTs. The distribution is fitted to a 
Gaussian plus an exponential and is normalized to unit area. The 0.2 PE threshold is clearly 
visible. 
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FIG. 10. Time response functions of the PMTs for two values of the predicted charge (a) [i 
= 0.0-1.0 PE and (b) [i = 5.0-6.0 PE for the isotropic + scattered Cerenkov light components in 
Michel electron events. The solid histograms show the data, and the points with error bars show 
the MC simulation. 
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FIG. 11. Corrected time distributions from "monoenergetic" data Michel electrons for (a) 
scintillation + scattered Cerenkov light and (b) direct Cerenkov light. The fit in (a) is to a 
convolution of a Gaussian with the sum of three exponentials and in (b) to a Gaussian. Both 
distributions are normalized to unit area. 
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FIG. 13. Cerenkov-to-scintillation density ratio p for cosmic-ray neutron events with elec- 
tron-equivalent energies between 60 and 200 MeV. 
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FIG. 14. Charge and time negative log-likelihoods for Michel electrons, as calculated for the 
entire event (top) and in the Cerenkov cone only (bottom) - analysis A. The solid histograms show 
the data, and the dashed histograms show the MC simulation. 
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FIG. 15. Charge and time negative log- likelihoods for "monoenergetic" Michel electrons, inside 
(top) and outside (bottom) the Cerenkov cone - analysis B. The solid histograms show the data, 
and the points with error bars show the MC simulation. 
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FIG. 16. Time difference distribution between the photons and the primary events in the 
initial DIF (beam on + off) data sample. The fit is to an exponential plus a constant. 
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FIG. 17. Cerenkov-to-scintillation density ratio p for all DIF data (beam on+off) for analysis 
B. Superimposed are the same distributions for cosmic-ray neutrons (dashed) and for the DIF-MC 
electron sample (dotted), normalized to the same area. 
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FIG. 18. Distributions of (a) the Fcer variable of analysis A and (b) the p variable of analysis 
B after all other selections have been applied. The solid histograms are DIF data (beam on + off), 
and the dashed histograms are for DIF-MC electrons, normalized to the same area. 



56 



UCRHEP-E191/LA-UR-97-1998 




FIG. 19. Charge and time negative log-likelihoods for all DIF data (beam on+off) for analysis 
A. Superimposed are the same distributions for the DIF-MC electron data (dashed), normalized 
to the same area. 
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FIG. 20. Charge and time negative log-likelihoods for all DIF data (beam on+off) for anal- 
ysis B. Superimposed are the same distributions for the DIF-MC electron data (dashed) and for 
cosmic-ray neutrons (dotted), normalized to the same area. 
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FIG. 21. Tank-veto timing correlations for all DIF data (beam on+off). The selected events 
are required to satisfy \t ve t \ > 50 ns in analysis A and \t veto \ > 70 ns in analysis B. 
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FIG. 22. Veto shield hit multiplicity distribution for the events in the final DIF sample of 
analysis B for the (a) beam-on, (b) beam-off and (c) beam-excess events. The solid histogram in 
(c) is the expected distribution from laser calibration events, which is due only to accidental hits 
in the veto system. 



60 



UCRHEP-E191/LA-UR-97-1998 



1 120 
> 

LU 

100 



80 
60 
40 
20 




180 



(a) 




100 200 300 400 500 600 700 800 

A (cm) 




700 800 
S (cm) 



FIG. 23. Distributions of (a) the A variable of analysis A and (b) the S variable of analysis B 
after all other selections have been applied. The solid histograms are DIF data (beam on + off), 
and the dashed histograms are for DIF-MC electrons, normalized to the same area. 
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FIG. 24. The distribution of event direction in the tank polar coordinate system denned by 
(cos 9 y , 4> xz ). The solid histograms show the beam-off DIF data, and the dashed histograms show 
the distributions for DIF-MC electrons (analysis A). 
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FIG. 25. The cos 6 V distribution, the cosine of the angle between the reconstructed event 
direction and that of the incident neutrino, for the final beam-excess DIF data events and that 
expected from the DIF-MC simulation (solid histogram). 
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FIG. 26. The reconstructed vertex to PMT surface distance distribution for the final 
beam-excess DIF events and that expected from the DIF-MC simulation (solid histogram). 
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FIG. 27. The x, y and z distributions for the final beam-excess DIF events and those expected 
from the DIF-MC simulation. 
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FIG. 28. Spatial distributions of the electron events of the final DIF sample in the x — y 
and y — z planes for (a) and (c) the 40 beam-on events and (b) and (d) the 175 beam-off events, 
respectively. The dotted contours outline the d > 35 cm fiducial volume. 
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FIG. 29. The energy distribution for the final beam-excess DIF events. The expectation for 
backgrounds (dotted histogram), the oscillation signal for large values of Am 2 (dashed histogram), 
and the sum of the two (solid histogram) are shown also. 
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FIG. 30. Photon to primary reconstructed distance for the old (dashed histogram) and new re- 
construction algorithms (solid histogram) for correlated photons. The average distance is decreased 
from 76 cm to 54 cm. The superimposed fit is to f(dr) = C dr 2 exp[— dr 2 /(2a 2 )]. 
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FIG. 31. Distribution of the number of photons with R y > 1 in the final DIF sample for the 
(a) beam-on, (b) beam-off and (c) beam-excess events (analysis B). 
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FIG. 32. The 95% confidence region for the DIF — > z/ e along with the favored regions for 
the LSND DAR measurement for &u —> v e . 
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